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The  phase  of  the  program  covered  by  this  report  was  devoted  to 
experiments  and  theoretical  modeling  on  discharge  pumping  of  the 
potassium-xenon  and  potassium-argon  excimer-dimer  systems. 

A time  dependent  discharge  model  was  developed  and  is  being  used  to 
predict  discharge  and  stimulated  emission  characteristics  o p 
transverse  discharge  through  potassium-argon  and  pota s sium-xenon 
mixtures.  Modeling  results  show  stable  discharge  operation  under 
conditions  which  yield  a practical  laser  gain  coefficient  on  both  the 
K-Xe  and  K-Ar  excimer  inds  and  on  the  K2  dimer  band,  j^ecif 
cally , gain  coefficients  on  the  order  of  10“  cm  are  P^e“\c_ted [JB 
the  K-Xe  or  K-Ar  excimer  bands  at  250  with  [Xe]  (or  [Ar])  ~ 10 
cm'3  at  a current  density  of  ~0.  1 A/cm3.  Similar  values  of  the 
gain  coefficient  on  the  K2  dimer  band  are  predicted  at  300  C IXel 
for  [Ar])  = 1 x 1020  cm-3,  and  at  a current  density  of  0.  n A/cm  . 

A discharge  apparatus  containing  transverse  electrodes  with  a 2 cm 
x 13.  5 cm  area  and  a 1 cm  gap  has  been  fabricated  which  is  capable 
of  operation  at  total  pressures  up  to  10  atm  and  a temperature  of 
350°C  Methods  for  obtaining  "aerosol-free"  mixtures  of  potassium 
and  argon  have  been  developed.  A diffuse  discharge  was  achieved 
with  a mixture  of  potassium  and  argon  at  100  psia  and  250  C.  D s- 
charge  electrical  and  spectral  measurements  were  carried  out  undei 
a variety  of  experimental  conditions  and  are  reported  herein. 
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I. 


INTRODUCTION 


This  is  the  third  semiannual  technical  report  on  the  ARPA/ONR 
program  to  develop  high  energy  alkali-rare  gas  excimer/dimer  lasers. 
The  primary  goals  of  this  program  are  to  identify  specific  operating 
conditions  for  achieving  laser  oscillation  on  the  A to  X transitions  of 
the  diatomic  alkali-xenon  excimer  molecules  (excimers)  and  of  the 
diatomic  homonuclear  alkali  molecule  (dimers),  and  to  demonstrate 
net  laser  gain  in  a laboratory  experiment  from  at  least  one  of  the 
alkali-xenon  systems. 

During  the  first  year  period  of  this  program  the  primary 
accomplishments  were  as  follows: 

1.  The  development  of  a theoretical  model  which  computes 
the  small- signal  gain  for  all  of  the  alkali-xenon  excimer/ 
dimer  systems  for  flashlamp  optical  pumping. 

2.  The  construction  of  a one  meter  long  high  temperature, 
high  pressure  optically  pumped  potassium  xenon  experi- 
mental cell. 

3.  Experimental  determination  of  the  absorption  coeffi- 
cient on  the  K-Xe  excimer  and  dimer  band  and  of 
the  stimulated  emission  coefficient  on  the  K-Xe  exci- 
mer band  all  of  which  were  consistent  with  the  modeling. 

4.  The  achievement  of  preliminary  observations  of  low 
pressure  discharge  pumped  fluorescence  spectra  of  the 
K-Ar  excimer  and  dimer  bands. 

Several  attempts  at  achieving  laser  oscillation  on  both  the  excimer  and 
dimer  bands  in  the  optically  pumped  K-Xe  tube  were  made  during  this 
first  year  period  and  during  the  early  part  of  the  present  reporting 
period.  Lasing  was  not  achieved  and  is  attributed  to  a combination  of 
factors  including  optical  loss  due  to  the  presence  of  aerosols  and  refrac- 
tive index  inhomogeneities,  a lack  of  adequate  flas^’amp  brightness, 
and  a possible  collisional  quenching  of  the  dimer  emission.  Time  did 
not  permit  any  of  these  factors  to  be  g~ven  adequate  study  in  the 
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optically  pumped  experiments  in  order  to  quantitatively  assess  their 
influence  on  the  net  laser  gain  before  the  shift  to  present  discharge 
pumping  experiments  was  implemented. 

Beginning  with  the  present  reporting  period  we  have  shifted 
from  optical  pumping  to  transverse  discharge  pumping  as  the  means  to 
demonstrate  efficient  scalable  laser  oscillation  of  an  alkali-rare  gas 
excimer/dimer  system.  Accomplishments  made  during  this  reporting 
period  include  (1)  the  development  of  a completely  time  dependent 
theoretical  model  for  discharge  pumping  of  the  K-Xe  excimer/dimer 
system,  (2)  the  construction  of  a transverse  discharge  tube  utilizing 
a 1 cm  x 2 cm  x 13.  5 cm  discharge  volume  and  a flashlamp  preionizer, 

(3)  the  apparent  elimination  of  aerosol  formation  in  this  tube,  and 

(4)  the  achievement  of  a diffuse  discharge  through  a potassium-argon 
mixture  at  ~200°C  and  a pressure  of  7 atmospheres.  Peak  currents  up 
to  300  A were  observed  at  a discharge  voltage  of  1 kV/cm  which  cor- 
responds  to  a discharge  power  loading  of  1 1 kW/cm  . The  current 
status  of  the  theoretical  and  experimental  activities  are  presented  in 
Sections  II  and  III,  respectively.  Section  IV  presents  the  anticipated 
work  plan  for  the  next  six  month  reporting  period  during  which  we 
anticipate  achieving  laser  oscillation  in  the  present  diffuse  discharge. 
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II. 


THEORETIC. 


-PROGRAM 


A,  Time  Dependent  Discharge  Model 

Both  the  optical  pumping  and  discharge  pumping  models  presented 
in  the  reports  prior  to  this  one  were  steady  state  models.  Steady  state 
modeling  is  justifiable  for  optical  pumping  since  there  the  relaxation 
times  of  all  optical  pumping  and  kinetics  processes  are  short  compared 
to  the  flashlamp  pulse  duration.  On  the  other  hand,  this  is  generally 
not  the  case  for  discharge  pumping.  In  particular,  such  discharges 
tend  to  be  unstable  and  it  is  not  a priori  clear  that  a steady  state  dis- 
charge mode  can  be  reached  at  any  time.  Thus,  the  original  steady 
state  discharge  model  discussed  in  the  previous  report  serves  only  as 
a preliminary  estimate  of  operating  conditions  for  achieving  laser 
threshold.  A full  time  dependent  discharge  model  is  presented  below. 

The  new  discharge  model  follows  the  time  evolution  in  the  dis- 
charge of  the  concentration  of  six  excited  and  ionized  alkali  states, 
the  electron  and  gas  temperatures , and  the  net  small-signal  gain  coeffi- 
cient at  two  wavelengths,  one  near  the  peak  of  the  excimer  gain  band, 
and  the  other  at  the  dimer  satellite  wavelength.  The  time  dependence 
of  the  circuit  current  and  voltages  is  also  included  with  a series  resist- 
ance, a series  inductance,  and  a storage  capacitor  as  the  modeled 
circuit  elements. 

The  plasma  and  gas  kinetic  processes  included  in  the  rate 
equations  governing  the  species  concentrations  are  listed  below.  The 
rate  constants,  cross-sections,  and  equilibrium  constants  for  these 
processes  can  be  found  in  the  computer  code  listing  includec  in 
Appendix  A,  and,  together  with  a source  reference,  in  Table  II- 1. 

An  R label  refers  to  a plasma  process  while  a G label  refers  to  a gas 
kinetic  process.  In  the  listing,  a B at  the  end  of  the  label  refers  to  the 
back  reaction  which  is  computed  by  detailed  balancing  using  the  appro- 
priate equilibrium  constant  or  Boltzmann  factor  reaction.  The  par- 
ticular case  presented  refers  to  a potassium- xenon  mixture. 
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Table  II- 1 . List  of  Cross-Sections,  Rate  Constants,  Equilibrium 
Constants  and  Source  References  used  for  Discharge 
Pumped  K-Xe/K^  Modeling 


Parameter 

Vaiue 

Reference 

Inelastic  electron  collision  cross- 
sections:  initial  slope  of  cross- 
sections  w.  r.  t.  energy  above 
threshold 

CPI 

CPI 

CP1I 

6 x ID'15  cm^/eV 
0.  5 x 10 ’ ^ cm^/eV 
5 x lO"1^  cm^/eV 

3 

3 

Rough  Estimate 

Recombination  coefficients: 
Three-body  RECP 
Dissociative  — (K  Xe+) 
RECPX 

Dissociative  — (K^) 

RECPD 

2.  28  x 10'26  x [K+]  x [Ne]  2 x Te('4'  39)  cm'3  sec'1 

3.  3 x 10‘5  x [K]  x [Xe+]  x [Ne)  x Te*'0- 67)  cm' 3 sec' 1 
3.  3 x 10* 5 x [K^]  x [Ne]  x Te('0-67)  cm' 3 sec'1 

5 (Cesium) 
Rough  Estimate 

Rough  Estimate 

Photoionization  Cross-Section: 
CPHI 

-19  2 

0.  5 x 10  17  cm 

4 

Electron  Elastic  Cross-Section: 
CEL 

5.  5 x 10'16  cm2 

3 

Heavy  Particle  Elastic  Cross-Section: 
CDIF 

1 x 10' 15  cm2 

Rough  Estimate 

Gas  Kinetic  Rates: 
GPD 

8 X 10'32  x [Xe]  X [K]  x [K*J  cm'3  sec*1 

6 

Equilibrium  Constants: 
KAEX 
KXDI 
KADI 

2.2  x 10"23  x exp  (0.074/Tg)  cm 
1 . 8 x 10'22  x exp  (0.  56/Tg)  cm3 
6.  5 x 10'23  x exp  (0.735/Tg)  cm3 

7 (Rubidium) 
6 
6 

Optical  Rates 

All  radiative  rates,  absorptive  and  stimulated  emission  cross-sections  and  their  dependences  on  temperature 
and  density  are  taken  from  Reference  2. 

— “ TIO-JA 
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The  plasma  processes  modeled  are: 


RPI: 

e"  + K 

— 

K+  + 2e~ 

RP1: 

e“  + K 

— 

* 

K + e 

RP1I: 

* 

e + K 

- 

K+  + 2e~ 

RECP: 

e"  + K+  + e" 

— 

K*  + e" 

RECPX: 

e"  + K Xe+ 

— 

K ' + Xe 

RECPD: 

e"  + *2  + 

- 

K*  + K 

The  gas  kinetic  processes  modeled  are: 

GPX: 

K + 2 Xe 

— 

K Xe  + Xe 

GPD: 

K + K + Xe 

— 

K2  + Xe 

GPX1: 

K*  + 2 Xe 

_ 

K Xe  ' + Xe 

* 

GPD1: 

K + K + Xe 

— 

K2  + Xe 

GPXI: 

K+  + 2 Xe 

— 

K XeV  + Xe 

GPDI: 

K+  + K + Xe 

K2+  + Xe 

The  radiative  processes  modeled  are: 

AT: 

j!c 

K 

- 

K + hv  (TRAPPED) 

AP: 

K 

— 

K + hv 

AP: 

KXe  ' 

- 

KXe  + hv 

AP: 

K2 

— 

K2  + hv 

Under  the  conditions 

we  consider 

([Xe]  « 

1 x 10  cm  j the 

gas  kinetic  rates  dominate  over  electron  collisional  rates  in  controlling 
the  population  of  the  molecular  states.  We  therefore  assume  thermal 
equilibrium  at  the  gas  temperature  to  be  maintained  between  the 


molecular  states  and  their  dissociation  products.  As  in  the  optical 

pumping  model  an  exception  is  the  ^ - state  of  K 2 for  which  radiative 
losses  are  included  since  they  can  compete  with  the  dissociation  rate. 

V he  average  electron  energy  is  computed  by  solving  the  electron 
energy  equation: 


J * Edisch 


= ^(INELASTIC  COLLISION  RATE  (Te))i  x E. 


+ (ELASTIC  COLLISION  RATE)  x 2 


(Me  \ 


(1) 


x Te 


where  J is  the  current  density,  Edigch  is  the  discharge  electric  field, 

E is  the  inelastic  energy  loss  for  the  ith  process  and  M ind  M are 
i me 

the  masses  of  the  M specie?  and  electron,  respectively.  At  present 
the  resonance  level  pumping  is  assumed  to  be  the  dominate  inelastic 
process  contributing  to  the  energy  equation.  For  this  reason  together 
with  the  fact  that  the  molecular  association  rates  which  pump  the  upper 
laser  levels  dominate  the  loss  of  the  resonance  level  population,  the 
predicted  efficiencies  for  converting  discharge  power  to  laser  power 
is  essentially  the  quantum  efficiency  of  the  transition  which  is  lasing. 

A Maxwellian  energy  distribution,  is  assumed  at  present  and  will 
eventually  be  replaced  by  an  energy  distribution  arrived  at  by  solving 
the  Boltzman  equation.  The  Maxwell  distribution  should  be  a fair 
approximation  because  of  the  dominance  of  two-step  alkali  ionization 
processes  over  the  single  step  alkali  ionization  under  the  considered 
operating  conditions.  Only  the  latter  process  requires  electrons  with 
energies  in  the  tail  of  the  distribution  where  a large  reduction  of  the 
number  of  electrons  from  a Maxwellian  distribution  is  known  to  occur. 
Excitation  and  ionization  of  the  rare  gas  is  neglected  entirely  for  this 
reason. 

The  electron  collisional  rate  constants  are  expressed,  as  in  the 
previous  report,  in  a form  which  results  from  assuming  a linear  rise 
of  the  cross-section  with  energy  above  threshold.  ^ 
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, b • 10  da 
R(Te)  - (6//rn  J/tt  de 


3/2 

Te  \ exp  (-e  /Te)  ( 1 + 


(x  4 zw) 


(2) 


where  € is  the  threshold  energy  for  the  excitation  or  ionization  process, 
da/de  £°  is  the  rat  of  rise  of  cross-section  with  energy  above  threshold 

and  m and  e is  the  electron  mass  and  change. 

6 Also,  in  solving  eq.  (1)  we  use  an  expression  for  the  drift 

velocity  obtained  from  Ref.  1 


= (2  T e/Mv  ) 


1/2 


(3) 


where  MXe  is  the  mass  of  the  xenon  atoms. 

The  gas  temperature  is  calculated  by  assuming 
discharge  power  goes  into  gas  heating  and  that  heat  is 
thermal  conduction  to  a wall  (assumed  cylindrical). 


that  all  of  the 
removed  only  by 


d TG  = j # E /(3/2  [M])  - (THERMAL  CONDUCTIVITY) 
<jt  disch'  , 

x (tjOTTEr)  x (Tg  ■ Twall) 

The  discharge  circuit  response  is  governed  by  the  familiar  circuit 
equations. 


dl 


disch 


= v.W-Uj-RM 


(5) 


vc  (t) 


charge 


1/c/ 


I dt 


(6) 


where  V ch  is  the  discharge  voltage,  Vc  is  the  capacitor  voltage, 
V is*  the  initial  charging  voltage,  I is  the  discharge  current, 

andR^L,  and  C are  the  series  resistance,  inductance,  and  storage 

capacitance,  respectively. 
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trSTIM(X)Ex  " (Ae/2c)  x X x REx  x dv  } Ex  (8) 

Di  Di  Di 

where  the  notation  follows  that  used  in  Ref.  2. 

As  shown  by  the  computer  code  listing,  the  coupled  rate  equations 
governing  all  of  the  above  processes  are  integrated  numerically  using 
a time  step  chosen  by  the  user  and  altered,  if  desired,  during  the 
course  of  the  computer  run.  The  computer  program  is  written  for  the 
PDP10  computer. 


B.  Results  and  Discussions 

The  results  of  three  computer  runs  of  the  above  model  are 
presented  in  Figs.  II- 1 through  II-3.  The  series  inductance,  L,  and 
resistance,  R,  are  assumed  equal  to  zero  in  these  runs. 

The  two  results  presented  in  Fig.  II- 1 are  for  an  avalanche 
discharge  through  a 10  atm  potassium-xenon  mixture  at  two  different 
temperatures.  As  in  the  optical  pumping  case,  the  lower  temperature 
result  shows  most  of  the  gain  occurring  on  the  excimer  band  (8500  A), 
while  the  higher  temperature  case  shows  the  highest  gain  occurring  on 
the  dimer  hand.  In  both  cases  the  model  predicts  that  a lairly  high 
E/N  value  of  ~10”  V/(cm-[K]  cm  ) is  required  to  give  practical  gain 
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coefficients  of  ~ 1%/cm.  The  much  lower  E/N  values  in  the  steady 
state  discharge  modeling  results  presented  in  the  last  semiannual 
report  were  due  primarily  to  the  erroneous  neglect  of  dissociative 
recombination  channels  for  the  alkali  ions  and  to  a lesser  extent  to  a 
poor  anproximation  used  therein  for  the  drift  velocity. 

The  results  in  Fig.  II- 1 are  for  an  infinite  storage  capacitance, 
and  illustrate  an  important  feature  of  the  discharge  kinetics.  Unlike 
the  other  popular  class  of  discharge  pumped  excimer  systems,  the 
rare  gas  monohalides,  the  alkaLi  rare  gas  discharge  appears  in  the 
m;deling,  to  be  stable  against  ionization  instabilities.  This  is  due  to 
two  causes:  first,  to  the  fact  that  the  discharge  is  recombination  stabil- 
ized rather  than  attachment  stabilized;  and  second,  to  the  fact  ground 
state  potassium  population  is  depleted  during  the  discharge  resulting 
in  a saturation  of  the  ionization. 

Results  for  a uv-sustained  K-Xe  discharge  are  presented  in 
Fig.  II-2.  Here  a comparable  dimer  gain  coefficient  is  achieved  at 
half  the  discharge  E/N  value  required  dor  the  300°  avalanche  discharge 
case.  Very  high  flashlamp  brightness  is  required  here  and  results  in 
an  extremely  low  overall  efficiency.  Since  the  higher  efficiency  avalanche 
discharge  is  predicted  to  be  stable,  there  is  no  apparent  advantage  in 
operating  the  discharge  in  this  mode. 

Finally,  in  Fig.  II - 3 we  present  modeling  results  for  a K-Ar 
discharge  with  a finite  storage  capacitance  at  conditions  near  those 
achieved  in  the  current  experiments.  The  potassium  density  is  seen  to 
be  too  low  to  yield  a practical  gain  coefficient  in  this  case.  We  also 
note  that  the  model  results  show  approximately  an  order  of  magnitude 
lower  peak  current  density  than  was  measured  in  the  experiment  at 
this  voltage  and  capacitance.  Tentatively,  we  attribute  this  to  too 
high  a value  estimated  for  the  dissociative  recombination  rates  used 
in  the  model.  We  plan  to  adjust  this  rate  to  match  the  measured 
discharge  characteristics  during  the  forthcoming  phase  of  experimental 
measurements. 
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RELATIVE  UNITS 


5197-2 


[Xe]  - 1 x 10^°  cm-3 
C * <*> 
v.  - 0 


lOkV/cm 


Edisch-40kV/cm 


GAIN  (10,400  A)  x 10  3 cm  ’) 
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Fig.  II- 1.  Theoretical  modeling  results  for  avalanche  sustained  dis 
charge  pumping  of  the  K-Xe/K^  system.  (a)Lowtemper 
ature  results,  (b)  High  temperature  results. 
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Fig.  II-2.  Theoretical  modeling  results  for  uv  sustained  dis- 
charge pumping  of  the  K-Xe/K^  system. 
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Fig.  II- 3.  Theoretical  modeling  results  for  avalanche  sustained  dis- 
charge pumping  of  K-Ar  system  under  conditions  achieved 
in  recent  experiments. 
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III. 


EXPERIMENTAL  PROGRAM 


The  primary  emphasis  of  the  experimental  part  of  the  program 
during  this  period  has  been  directed  toward  the  achievement  of  a diffuse 
discharge  in  potassium-argon  mixtures  at  high  pressure.  The  dis- 
charge dimensions  and  configuration  were  chosen  in  order  to  obtain 
experimental  results  which  are  directly  applicable  to  scaling  consi- 
derations for  high  power  laser  applications.  Thus  the  discharge  ap- 
paratus was  designed  around  a transverse  electrode  arrangement  with 
the  width  to  gap  ratio  equal  to  two  and  the  length  to  gap  ratio  equal  to 
ten. 

In  addition  to  design  and  fabrication  of  the  electrode  apparatus, 
considerable  attention  was  devoted  to  solution  of  the  "aerosol"  problem 
encountered  in  our  previous  work.  Experiments  were  conducted  in  a 
carefully  controlled  environment  under  a variety  of  experimental  con- 
ditions which  have  demonstrated  that  "aerosol-free"  mixtures  of  po- 
tassium and  argon  can  be  obtained  at  elevated  temperatures  and  pres- 
sures provided  temperature  gradients  are  minimized  and  the  potassium 
vapor  is  evolved  at  low  pressure.  Some  of  our  experiments  indicate 
that  the  former  requirement  is  more  important  than  the  latter. 

During  the  latter  part  of  this  reporting  period  we  have  achieved 
a diffuse  discharge  in  a potassium-argon  mixture  at  high  pressure 
(7  atm).  The  details  of  these  first  experiments  will  be  discussed  in 
the  following  section  but  it  is  noteworthy  that  this  discharge  has  very 
good  spatial  homogeneity,  exhibits  strong  K-Ar  excimer  emission  and 
that  no  form  of  auxiliary  preionization  was  required.  Potassium  atom 
densities  utilized  to  date  are  not  high  enough  to  provide  significant  op- 
tical gain;  however,  apparatus  modifications  are  presently  underway 
which  will  allow  the  achievement  of  equilibrium  potassium  vapor  densi- 
ties at  temperatures  up  to  350°C.  Recent  measurements  indicate  that 
inadvertent  heat  transfer  through  mounting  apparatus  has  restricted  the 
maximum  oven  temperature  to  250°C. 
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A. 


Discharge  Apparatus 


A major  concern  in  the  apparatus  design  was  the  electrode  size 
and  configuration.  An  electrode  size  of  2 cm  x 13.  5 cm  with  a gap 
spacing  of  1 cm  was  chosen  initially  because  it  is  small  enough  to  be 
convenient  but  also  large  enough  to  give  meaningful  results  in  terms  of 
scaling  the  system  for  high  power  laser  applications.  In  order  to  allow 
maximum  flexibility  in  choice  of  discharge  mode  for  optimum  laser 
operation  we  have  included  provisions  in  the  apparatus  for  a high  inten- 
sity flashlamp.  This  provides  the  option  of  discharge  operation  in  a uv 
preionized  mode,  a uv- sustained  mode  or  simply  as  an  over- voltaged 
avalanche  discharge  without  preconditioning.  A diagram  of  the  first 
version  of  this  system  is  shown  in  Fig.  Ill- 1 . A pyrex  pressure  vessel 
was  selected  to  provide  maximum  viewing  capability  of  the  discharge. 
This  design  proved  to  be  too  difficult  to  fabricate,  however,  so  an  al- 
ternate apparatus  utilizing  a stainless  steel  pressure  vessel  with  a 
pyrex  liner  for  electrical  isolation  was  designed  and  assembled.  This 
apparatus  is  illustrated  in  Fig.  Ill- 2 where  the  electrodes,  re-entrant 
Brewster  angle  windows,  and  a large  side  window  for  viewing  the  dis- 
charge along  its  entire  length  are  shown.  This  apparatus  met  the  high 
pressure  (10  atm)  and  discharge  viewing  requirements  but  problems 
developed  with  the  "O"  ring  seals  after  several  temperature  cycles 
between  25°C  and  350°C.  Chemical  reactions  between  several  of  the 
components  caused  difficulties,  and  thermal  gradients  contributed  to 
aerosol  formation.  This  system  was  revised  with  the  use  of  metal 
vacuum  flanges  and  elimination  of  the  pyrex  liner.  A photograph  of 
this  apparatus  is  shown  in  Fig.  Ill- 3.  Chemical  corrosion  problems 
were  virtually  eliminated  as  well  as  aerosol  formation  in  the  new  appa- 
ratus. A compromise  was  made  on  discharge  viewing  length  because 
of  the  maximum  window  size  available  with  a metal  seal;  additional 
windows  can  be  added  for  better  discharge  viewing  if  required  in  later 
experiments.  This  is  the  device  with  which  we  have  observed  spatially 
uniform  discharges  in  a potassium-argon  mixture  at  7 atm.  The  major 
problem  encountered  with  this  design  has  been  a limitation  on  maximum 
temperature  of  ~ 250°C  caused  by  heat  loading  at  the  mounting  fixtures. 
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Fig.  Ill  - 1 . Alkali-rare  gas  discharge  tube,  design  - 
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Fig.  Ill  - 3 . Photograph  of  alkali-rare  gas  discharge  appara- 
tus, design  3. 
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Modifications  are  presently  underway  to  correct  this  design  limitation. 

The  electrode  profile  v'hich  we  have  used  to  date  is  planar  with 
the  sides  and  ends  contoured  to  the  curvature  of  an  ellipse.  Ihe  elec- 
trode is  flat  over  a region  ~ 1.  3 times  the  discharge  gap  and  then  joins 
smoothly  to  an  ellipse.  In  previous  work  we  have  found  the  electric 
field  between  electrodes  of  this  type  to  be  uniform  to  better  than  1% 
as  determined  by  electrostatic  potential  plots  using  resistance  paper. 
The  present  electrode  profile  is  expected  to  be  satisfactory  as  judged 
by  our  first  results  and  is  advantageous  for  scaling  to  larger  sizes 
because  of  the  simplicity  of  its  shape.  The  major  difficvlty  with  other 
electrode  shapes  such  as  the  Rogowskii  profile  is  the  extreme  width 
which  must  be  used  if  the  central  portion  of  the  electrode  is  made  very 
planar  to  avoid  constriction  of  the  discharge  width. 

B.  Aerosol  Studies  and  Optical  Measurements  of  Potassium 

In  some  of  our  previous  work  with  small-scale  apparatus 

3 

(75  cm  ) we  have  observed  significant  optical  scattering  losses  in 
potassium- rare  gas  mixtures  at  high  temperatures  (300  C)  and  pres- 
sures (10  atm).  This  was  also  found  to  be  a problem  with  the  large 
volume  apparatus  shown  in  Fig.  Ill- 2,  Quantitative  optical  transmis- 
sion measurements  and  qualitative  side  scattering  observations  were 
conducted  during  this  reporting  period  with  apparatus  of  the  type  shown 
schematically  in  Fig.  III-4.  Results  of  transmission  measurements 
with  a broad-band  tungsten  light  source  are  illustrated  in  Fig.  III-5. 

It  is  apparent  from  the  data  in  Fig,  III- 5 that  a large  optical  loss  occurs 
over  a wide  wavelength  region  ( = 400  A)  under  certain  experimental 
conditions.  In  addition  to  this  transmission  reduction  with  time,  strong 
side  scattering  from  a He-Ne  laser  beam  was  also  observed.  This 
experiment  was  conducted  with  the  apparatus  shown  in  Fig.  III-2. 
Potassium  was  vaporized  from  a 0.  6 cm  x 15  cm  stainless  steel  wick 
situated  near  the  bottom  of  the  housing  with  an  argon  pressure  of  1 atm. 
Following  these  observations  of  aerosol  formation  in  the  discharge 
apparatus,  diagnostic  experiments  were  conducted  in  sealed  off  pyrex 
vessels  situated  in  a uniform  temperature  muffle  furnace  as  shown  in 
Fig.  HI-6. 
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Fig.  III-5.  Optical  tiansmission  data  illustrating  aerosol  for- 
mation in  potassium-argon  mixture. 
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Optical  transmission  measurements  with  the  broad-band  light 
source  over  400  A bandwidths  centered  at  7700  A and  6475  A at  different 
pressures  of  argon  were  carried  out  along  with  transmission  and  scat- 
tering measurements  with  a He-Ne  laser  beam.  The  cell  was  filled 
with  potassium  and  argon  on  a separate  processing  apparatus,  sealed 
off  and  located  in  the  furnace  for  optical  measurements.  Quantitative 
correlations  with  equilibrium  atomic  potassium  and  potassium  dimer 
densities  were  observed  over  the  temperature  range  of  25  C to  350  C 
and  argon  pressures  from  0 to  2 atm;  experiments  were  not  conducted 
at  higher  argon  pressures  in  this  apparatus.  A comparison  between 
our  measured  values  of  potassium  dimer  absorption  coefficient  and 
those  of  Lapp  and  Harris8  is  shown  in  Fig.  III-7.  Using  their  value  of 
20  A^  for  the  absorption  cross  section  of  at  6465  A and  our  mea- 
sured absorption  coefficient  of  0.098  cm  1 at  325  C,  we  obtain  a potas- 
s ium  dimer  density  of  4.  9 x 10°  cm"3  which  compares  very  favorably 
with  the  density  of  5.0  x 1013  cm'3  calculated  from  equilibrium  vapor 
pressure  data^.  It  is  more  difficult  to  deduce  atomic  densities  of  po- 
tassium from  optical  absorption  measurements  at  high  temperatures 
and  pressures  because  of  the  complex  optical  broadening  mechanisms 
involved;  a well  documented  curve  of  growth  for  the  integrated  absorp- 
tion under  these  conditions  has  not  been  established.  However,  since 
we  have  obtained  good  agreement  between  equilibrium  values  of 
potassium  dimer  densities  and  our  measurements,  it  is  reasonable  to 
infer  that  equilibrium  between  atomic  and  dimeric  potassium  is 
attained  in  our  apparatus.  Thus  we  will  obtain  atomic  potassium 
densities  from  measurements  of  dimeric  potassium  under  conditions 
where  the  dimer  spectrum  is  pressure  broadened  into  a continuum. 

No  aerosol  formation  was  observed  during  the  experiments 
conducted  in  this  apparatus.  This  is  illustrated  by  the  optical  trans- 
mission data  shown  in  Fig.  Ill- 8 which  was  obtained  under  the  care- 
fully controlled  conditions  described  above.  Note  the  dramatic  contrast 
between  this  data  and  that  shown  previously  in  Fig.  III-5.  We  conclude 
that  equilibrium  values  of  potassium  vapor  density  can  be  obtained  in 
argon  at  temperatures  up  to  350°C  and  pressures  up  to  2 atm  if  the 
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Fig.  III-7 . Potassium  dimer,  (a)  Optical  transmission,  and 
(b)  absorption  coefficient.  Solid  curve  is  from 
Feference  8,  dashed  curve  is  our  measurement. 
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Fig.  Ill  - 8 . Optical  transmission  through  potassium-argon  in  15  cm 

sealed  off  cell. 
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temperature  of  the  apparatus  is  kept  uniform.  Ihis  is  considerably 
more  difficult  in  the  discharge  system  without  incorporating  elaborate 
heating  equipment  because  of  the  uneven  heat  capacity  of  the  internal 
parts  of  the  discharge  chamber.  We  infer  that  potassium  aerosols 
form  when  macroscopic  size  volumes  of  saturated  vapor  densities  come 
in  contact  with  cold  spots  inside  the  discharge  chamber.  In  the  region 
around  these  cold  spots  the  saturated  vapor  density  is  transformed  to  a 
super- saturated  vapor  through  heat  transfer  to  the  colder  region.  The 
super- saturated  vapor  is  subsequently  converted  to  an  equilibrium  va- 
por by  condensation  into  macroscopic  particles  of  liquid  potassium. 
High  pressures  of  inert  gas  in  a large  diameter  chamber  contribute  to 
the  aerosol  formation  by  increasing  the  diffusion  time  of  the  alkali  va- 
por and  thereby  preventing  it  from  maintaining  a dynamic  equilibrium 

with  the  coldest  regions  in  the  chamber. 

With  this  hypothesis  and  our  experience  with  the  sealed- off 
cells  in  the  uniform  oven,  we  postulated  that  aerosols  could  be  elimin 
ated  in  the  large  discharge  apparatus  by  slowly  evolving  the  potassium 
vapor  in  the  absence  of  inert  gas  and  then  mixing  pre-heated  gas  with 
the  alkali  vapor  after  equilibrium  vapor  densities  of  tne  pure  metal 
vapor  were  obtained.  Experiments  were  performed  following  this  pro- 
cedure by  pre-heating  argon  to  a temperature  50  degrees  higher  than 
the  temperature  of  the  discharge  chamber.  Uniform  mixing  was 
achieved  by  admitting  the  gas  into  the  potassium  vapor  by  means  of  a 
perforated  stainless  steel  sting.  The  results  of  this  set  of  experiments 
are  summarized  by  the  optical  transmission  data  in  Fig.  Ill- 9 taken 
with  the  apparatus  as  shown  schematically  in  Fig.  Ill- 4.  No  aerosol 
formation  occurred  as  indicated  by  this  data  and  by  the  complete  ab- 
sence of  scattering  from  the  He-Ne  laser  beam.  Thus,  we  have 
arrived  at  an  experimental  procedure  for  obtaining  homogeneous  mix- 
tures of  alkali  vapors  and  rare  g .ses  at  high  temperatures  (300°C)  and 
pressures  (7  atm).  These  values  of  temperature  and  pressure  are  by 
no  means  expected  to  be  upper  limits  out  rather  reflect  the  range  over 
which  carefully  controlled  experiments  have  been  conducted  in  appara- 
tus of  the  type  necessary  for  alkali-rare  gas  exdmer  laser  systems. 
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Fig.  Ill- 9 - 


Optical  transmission  through  the  discharge  ap 
atus  indicating  no  aerosol  formation.  (T  = 300 


ar- 

C). 


31 


C.  Electrical  Discharge  Experiments 

After  establishing  an  aerosol-free  mixture  of  potassium  vapor 

and  argon  at  7 atm  pressure,  electrical  discharge  measurements  were 

19-3 

initiated.  Using  this  argon  density  (9x10  cm  ) and  an  indicated 
oven  temperature  of  300°C,  a variety  of  experiments  were  conducted 
to  evaluate  both  the  electrical  and  optical  characteristics  of  the  dis- 
charge. Qualitative  visual  observations  of  the  discharge  through  the 
1. 5 cm  side  viewing  window  and  through  the  end  windows  indicated 
that  we  had  achieved  a very  uniform  diffuse  discharge  over  the  entire 
electrode  area  without  the  use  of  any  discharge  preconditioning. 

This  is  partially  illustrated  by  photographs  of  the  discharge 
through  the  side  window  shown  in  Fig.  Ill- 10.  The  field  of  view  of  the 
camera  used  here  is  considerably  smaller  than  the  full  discharge  re- 
gion and  also  smaller  than  a person's  eye  located  adjacent  to  the  side 
window.  Two  exposure  levels  are  shown  to  verify  the  lack  of  small 
scale  striations  in  the  discharge.  As  far  as  can  be  ascertained,  the 
discharge  is  very  uniform  under  the  operating  conditions  of  temperature 
and  pressure  utilized  here  (300°C,  7 atm  Argon). 

Voltage  and  current  characteristics  of  this  pulsed  discharge 
operating  without  preionization  are  shown  in  Figs.  Ill- 11  and  III-12. 

A negative  voltage  is  applied  across  the  electrodes  by  .neans  of  a thy- 
ratron  switching  tube  and  a capacitor  previously  charged  to  a particu- 
lar voltage  level.  The  current  pulse  develops  after  a delay  time  of  a 
few  microseconds  and  has  a time  duration  which  is  dictated  by  the  self- 
inductance and  capacitance  of  the  circuit  rather  than  the  plasma  im- 
pedance. As  the  applied  voltage  is  increased,  breakdown  occurs  at 
earlier  times  and  the  peak  current  increases  in  proportion  to  the  charge 
initially  stored  on  the  capacitor.  This  variation  in  time  between  appli- 
cation of  the  voltage  and  development  of  the  current  is  not  visible  in 
Figs.  HI- 11  and  III-  12  because  the  complete  time  interval  of  the  experi- 
ment can  not  always  be  displayed  on  the  time  scale  chosen  to  fully 

resolve  the  current  pulses.  Considerable  power  loading  into  the  dis- 

3 

charge  has  been  achieved  (~10  kW/cm  ) without  development  of  an  arc. 
When  the  temperature  is  decreased  by  25  C from  the  indicated 
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Fig.  II I - 10.  Photograph  of  potassium-argon  discharge 
as  viewed  through  the  1.  5 cm  side  win- 
dow. Film  exposure  of  (b)  is  5 times 
less  than  (a). 
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V0  = 8 kV 
C *=  2.7  nF 


Fig.  Ill— 11.  Current-voltage  waveform 

characteristics  of  potassium- 
argon  transverse  discharge. 
Upper  traces  — voltage, 
lower  traces  — current. 
Sensitivities:  (a)  5 kV/cm, 

20  A/cm;  (b)  5 kV/cm, 

50  A/ cm. 
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Vq  = 5 kV 
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V0  = 8 kV 
C =0.01  gF 


Fig.  Ill- 12.  Current-voltage  character- 
istics of  potassium-argon 
discharge.  Upper  traces  — 
voltage,  lower  traces  — 
current.  Sensitivities: 

(a)  5 kV/cm,  50  A/cm; 

(b)  5 kV/cm,  200  A/cm. 


55 


temperature  of  300°C,  arcing  does  occur  approximately  50%  of  the 
time  and  arcs  were  observed  occasionally  at  voltages  near  9 kV.  How- 
ever, the  latter  case  may  be  simply  related  to  an  internal  breakdown 
through  the  high  voltage  feedthrough  to  ground  rather  than  across  the 
electrodes.  Minor  changes  in  the  apparatus  are  presently  being  made 
in  order  to  prevent  this  possibility  in  future  work.  A few  experiments 
were  also  carried  out  using  a quartz  flashlamp  for  pre-iomzation  of 
the  discharge.  The  results  have  not  been  examined  in  detail  but  we 
note  here  that  the  time  interval  between  discharge  current  onset  and 
voltage  application  decreases  markedly  (several  microseconds)  with 
increasing  flashlamp  intensity. 

In  addition  to  characterization  of  the  electrical  features  of  the 
rh" charge,  we  have  also  conducted  spectroscopic  investigations  of  the 
cptical  emission  at  wavelengths  from  7400  A to  1.1pm.  An  emission 
jet  trum  in  the  wavelength  region  near  the  potassium-argon  excimer 
hand  is  shown  in  Fig.  III-13(a).  Self-absorption  by  ground  state  po- 
tassium in  the  region  between  the  electrodes  and  Brewster  angle  win- 
dows  is  responsible  for  the  two  negative  peaks  at  7665  A and  7669  A. 
The  length  of  this  optical  loss  region  will  be  reduced  from  the  present 
5 cm  to  ~ 1 cm  during  the  device  modifications  that  are  currently  in 
progress.  Another  spectrum  is  shown  in  Fig.  III-13(b)  with  an  inten- 
sity scale  2.5  times  more  sentitive  than  in  Fig.  III-13(a)  which  demon- 
strates the  large  bandwidth  and  spectral  shape  of  the  excimer  emission 
The  emission  peak  at  808  5 & has  not  been  specifically  identified,  but 
may  be  a result  of  impurities  in  the  system.  No  emission  could  be 
detected  at  wavelengths  greater  than  8700  K although  strong  potassium 
dimer  emission  is  expected  at  1 . 04  pm.  After  a thorough  search  for 
the  dimer  emission,  we  dismantled  the  discharge  apparatus  and  found 
that  potassium  metal  had  condensed  near  one  of  the  system  mounting 
posts.  This  suggests  that  the  coldest  position  in  the  apparatus  was 
at  this  location  and  was  significantly  lower  in  temperature  due  to  heat 
transfer  through  the  mounts  than  indicated  by  the  probe  thermocouple 
located  in  a quartz  tube  near  the  center  of  the  electrode  regions.  Ten 
other  thermocouples  were  located  throughout  the  apparatus  and  these 
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Fig.  Ill  - 1 3 . Potassium-argon  excimer 

emission.  Sensitivity  of  (b)  is 
2.  5 times  that  of  (ah  Wave- 
length scale  — 100  A/cm. 
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positions  were  always  mai  itained  at  higher  temperatures  than  the  probe 
thermocouple  in  the  quartz  tube.  Moreover,  after  these  discharge 
measurements  were  made,  we  carefully  monitored  the  optical  trans- 
mission of  a He-Ne  laser  beam  through  the  discharge  apparatus.  This 
data  is  presented  in  Fig.  Ill- 14.  In  addition  to  a residual  optical  loss 
at  the  end  windows,  the  onset  of  absorption  occurs  at  a measured 
temperature  which  is  50°  higher  than  the  calculated  temperature  for 
absorption  onset  based  on  equilibrium  vapor  pressure  data.  The  actual 
or  effective  gas  temperature  is  clearly  considerably  lower  than  indi- 
cated by  the  thermocouple  in  the  quartz  probe  tube.  Thus  we  conclude 
at  this  time  that  the  discharge  measurements  to  date  in  this  apparatus 
have  been  made  at  temperatures  which  may  be  bO°C  or  more  lower 
than  measured  directly.  This  does  not  reduce  the  significance  of  the 
results  obtained  but  mandates  equipment  modifications  before  potassium 
vapor  densities  high  enough  to  achieve  optical  gain  can  be  obtaine' 

We  are  now  in  a position  to  accurately  monitu.  atomic  potassium 
through  optical  measuj  - mcnts  of  potassium  dimer  concentrations  since 
we  have  experimentally  demonstrated  a correlation  with  our  measure- 
ments and  equilibrium  vapor  pressure  data  of  both  these  species  (see 
Fig.  IT.I-7  and  the  discussion  in  Section  III-B). 

An  additional  experimental  observation  on  the  discharge  emis- 
sion is  shown  in  Fig.  Ill- 15  where  it  is  seen  that  the  excimer  emission 

O 

in  the  7700  A region  closely  follows  the  temporal  behavior  of  the  dis- 
charge current  pulse.  In  Fig.  Ill- 16  is  shown  a comparison  between 
the  potassium-argon  excimer  emission  and  absorption  which  illustrates, 
along  with  Fig.  Ill- 13(b),  the  extension  of  the  long  wavelength  excimer 
emission  beyond  the  absorption  region.  It  is  in  this  spectral  region 
where  the  excimer  emission  is  strong  compared  to  absorption  and  con- 
sequently where  optical  gain  can  be  obtained. 
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TRANSMISSION. 


5107-18 


Fig.  Ill- 14.  Potassium  dimer  concentration  and  optical  transmission  at 
6328  K. 
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Fig.  Ill-  1 5. 

Temporal  profiles  of 
potassium-argon  excimer 
emission  (upper  trace)  and 
discharge  current  (lower 
trace). 
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7699  A 


TRANSMISSION 


EMISSION 


Fig.  Ill- 16.  Superposition  of  potassium-argon  excimer 
optical  emission  and  transmission  spectra. 
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IV. 


FUTURE  WORK  PLAN 


All  theoretical  and  experimental  efforts  during  the  next  six-month 
period  will  be  directed  toward  achieving  laser  oscillation  on  either  the 
K-Xe  or  K-Ar  Z-X  excimer  transition  or  the  A-X  d'-ner  transition 
from  the  present  transverse  discharge  apparatus. 

Modifications  to  the  experimental  discharge  apparatus  which 

16  - 3 

will  allow  the  achievement  of  potassium  densities  up  to  10  cm  are 
presently  in  progress.  Additionally,  a new  set  of  electrodes  having 
a length  of  25  cm  is  being  fabricated  which  represents  an  increase  in 
optical  gain  path  of  11.  5 cm  over  the  previous  electrodes. 

Theoretical  modeling  will  concentrate  on  guiding  the  experiments 
into  operating  regimes  which  appear  to  provide  the  highest  laser  gain 
under  conditions  generally  being  achieved  in  the  laboratory.  Up-grading 
of  the  modeling  will  be  carried  out  continuously  through  the  incorpora- 
tion of  cross-section  and  rate  constant  data  being  obtained  or  inferred 
through  measurements  here  and  elsewhere.  Also,  the  Maxwellian  elec- 
tron distribution  used  at  present  will  either  be  quantitatively  justified 
or  replaced  by  either  a solution  to  the  Boltzmann  equation  or  a two- 
group  Maxwellian  distribution.  Stimulated  emission  processes  will 
be  added  to  the  model  when  appropriate  for  interfacing  with  the 
experiments. 
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APPENDIX  A 


C KXECP  - PULSED  K-XE  LASER 

REAL  MASE,MASXE,MASM,M1 ,M2,MI ,K,K3X,K5X,K4I ,NEO,NE 

REAL  LAMDE , LAMDD , KAEX , KADI , KXEX , KXD I 

REAL  NEP 

REAL  J 

REAL  LENGTH 

REAL  KPLUS.KNI ,KEL 

REAL  L 

EXTERNAL  FUNC 

CO  MMON/STOR/XE , PO  f CP  I ,CXI ,NE,NEPtL,A,R,DT,E, 

+ EPI ,EXI , CEL, MASE.MASXE.V, LENGTH, PI ,CPI I,CPI,EPI 
COMMON/STE/TE 

C CONSTANT  DATA  (DATA  IS  FOR  K-XE) 

MASE=0.5E6l  E=I.6F-19»  MASXE= I 3 I * 1 840*. 5E6 
R=0 

C=2 • 7E-9 
D=  I 
A=27 

LAMDE=.85E-4 
LAMDD= I . 04E-4 
AP=3.69E7 
REX=3 . 7E-8 
RDI=5.4E-8 
RDIE=4. 4E-8 
DRDFE=. 53E-8/( 579 ) 

DRDFD= I E-8/( .04E4) 

DRDFDE=.4E-8/( I . 388E3 ) 

VEA=-.079 

VEX=.062 

VDA=-.66 

VDX=-. 262 

VDXE=-.560 

EPI  = I .61 

EPI=4. 34 

EPX I = I . 46 

EPX  1 = 4.  19 

E PD  I = 1 .46 

EPDI=4. 19 

CP  I =6E- 1 5 

CPI=.5E-I 6 

CPI  I = 5E - 1 6 

CPH I=.05E- 1 8 

CEL=5 .5E- 1 6 

CDI F= 1E-I 5 

SA=8. 49E4?SB=7. 18 

LENGTH= I 

K= I 60E I 9 
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AX2=4E7 ; AP=3.69E7  tAXEDI =4E7 
CSTME=AP/2*LAMDE**2*REX**2*DRDFE 
CSTMD«AP/2*LAMDD**2*RDI**2*DRDFD 
CSTMDE=AP/2*LAMDE**2*RD1'E**2*DRDFDE 

C INITIAL  CONDITIONS 
WRITE! 5 , 50) 

50  F()RMAT(2X,2HXE,2X,2HTG,2X, I HW, 2X, I HV, 2X , IHC,2X,2HDT, 2X,6HIPRINT, 
+2X,4HTMAX, 2X,5HDELAY) 

ACCEPT*, XE, TO, W,V,C,DT, I PRINT,TMAX, DELAY 
P0=2.7E I6*( 273/(273+TG) ) *1 0**(SB-.052*SA/(TG+2 73 ) ) 

TG= (TG+273)*! 1/1 . 16E4) 

TGO=TG 

F=W/( I . 6E-I 9*EP I ) *( 1.24/EPI-.2) 

KXD I = 1 .8E-22*EXP( . 56/TG) 

KXE  X=2 . 2E-23*EXP ( -. 05/TG ) 

PX=KXEX*PO*XE 
PD=KXDI *P0**2 
PX=KXEX*P0*XE 
PI-F*CPHI*PO*DELAY 
PI=10.*PI 

AT=AP*I . 6/(3. 3 E- I 5*P0*D) 

IF(AT-AP)  13,17,17 

17  AT=  AP 

18  CONTINUE 

WRITE(5,  82)  PO, PD, PI 

82  FORMAT(  2X,  /P0=/ , E 1 0 .3, 5X,  /PDss/,  E I0.3,5X,/PI  = /,EI0.3) 

WRITE (5,  80) 

80  FORMAT( 4X, 4HTIME ,6X , 5HVDISC , 5X , 1 HJ,9X,2HNE ,8X , 2HPI , 

+8X,3HPXI ,7X,3HPDI ) 

WRITE (5,  85) 

85  F0RMAT(7X, 1HP,9X,2HPI ,8X,5HGAINE,5X,5H0AIND,5X,5HBETAE,5X,5HBETAD, 
+5X, 2HTE ) 

C TIME  STEP 

20  ISTEP=ISTEP+1 
T=T+DT 

IF ( T-TMAX)  21 ,21 ,22 

21  CONTINUE 

C CIRCUIT  RESPONSE 
P=PO-(PI+PI) 

NE=PI 

V=V-DT* J*A/C 

IERR=0; IND= I ;EPS=0.0l 

TE«SOLn(IND,FUNC,. I ,5., EPS, I ERR) 

IF(IERR)  1,2,1 

1 WRITE(5,*)  I ERR; GO  TO  30 
STOP 

2 CONTINUE 
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C HATE  CONSTANTS  - ELCTRONIC 

V D=  SQ  RT  ( 2 *TE  /M  AS  XE ) *3E I 0 
J«NE*E*VD 

VDISC=V-R*J*A-L*A*E*VD*(NE-NEP)/DT 


NEP=NE 

RC=6.38E7*(SQRT(2*TE ))**3 

RPI =RC*EXP( -EPI /TE )★( I +EP1/(2*TE) )*(NE*CPI )*P 
RPI B=RP I /P*EXP(EPI /TE)*P I *. 33 
RPI=RFN(EPI )*(NE*CPI )*P 
RPI I*RrN(EPI-EPI )*(NE*CP I I )*PI 
RECP=PI *2. 28E-26*NE*TE**( -4.39) +NE 
RECPX=3.3E-5*TE**(-.67)*PXI*NE 
RECPD=3.3E-5*TE**(-.67)*PDI*NE 
C RATE  CONSTANTS  - GAS  KINETIC 
SQRTG=SORT(TG) 

KAEX=2.2E-23*EXP( .074/TG) 

KADI =6. 5E-23*EXP( . 735/TG) 
KXEX=2.2E-23*EXP(-.05/TG) 

KXDI = I .8E-22*EXP( . 56/TG) 

GPX=XE*8E-32*XE*P I 
GPXB=XE*8E-3 2/KAEX 
GPD=P*( XE*8E-30*P I ) 

GPDB=(XE*3E-30)/KADI 
GPXI=(XE*8E-32*XE)*PI 
GPXIB=(XE*8E-32*PXI )/KAHX 
GPDI=(XE*8E-30*P)*PI 
GPDIB=< XE*8E-30*PDI )/KADI 


C DIFFUSION  RATES 

KP! US=2E3*( 3E I 6/( XE ) )*SGRT( . 03/TG) 

KEL=VD/ ( VDI SC/LENGTH ) 

DXE=SQRT(2*TG/MASXE )*3E I 0/( 3*( XE)*CDIF)*( 4. 8/D) **2 
DAPOS=KPLUS*TE*( 4.8/D) **2 
DP=DXE* I .8 

WA<lt"p?aN+l)T*(RPI+HECP+Hb'CPX+HECPD-RPlli-AP*(PXI+PDI  )-DP*PI  ) 

IF  (P1/P-3*EXP(-EPI/TE) ) 8,9,9 
9 p l=P*3*EXP(-EPI/TE) 

8 CONTINUE 
PX=KXEX*P*XE 
PD=KXDI*P**2 

PI=PI +DT*( RPI+RP 1 1 +RP2 I+GPN*P+P*CPHI*F-RECP-hECPX-RECPD-DAPOS*P I 
PX I =GPX/( GPXB+AP ) 

PXI=KAHX*PI*XE 
PD 1 -GPD/(GPDB+AP) 

PDI=KAD!*P*PI 

C uAS  ™^™^(^H^VL)^VDISC/Ll:NGTH_(4#8/D)**2*K*(TG-TG0))*2/(3*XE) 


) 
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C GAIN 

BETAE«CSTME*2*P*XE*EXP(-VEX/TG)/3EI0+CSTMDE*P**2*EXP(-VDXE/ 

+TG) ★. 25/3E 1 0 

BETAD«CSTMD*P**2*EXP(-VDX/TG)*.25/3EI0 

GAINE*CSTME*(PXI*EXP(-VEA/TG)/(KAEX*1 .5 ) -2*P*XE*EXP( -VEX/TG ) )/3E 1 0 
+-C3TMDE*(P**2*EXP(-VDXE/TG)*.25)/3EI0 
GAIND«CSTMD*<PDI*EXP(-VDA/TG)/<KADI*12)-P**2*EXP(-VDX/TG)*.25) 
+/3E10 
C OUTPUT 

IF C ISTEP/IPR INT-FLOATC I STEP) /FLOAT( IPRINT) ) 31 ,30,31 

30  CONTINUE 
WRITE (5,*)  TG 

WRITE (5,90)  T,  VDISC, J,NE, PI ,PX I , PD  1 
90  E0RMAT( 7CE10.2) ) 

WRITE ( 5 , 95)  P,P 1 ,GAINE ,GAIND, BETAE,BETAD,TE 
95  FORMAT ( 3X, 7CE10.2 ) ) 

31  CONTINUE 
GO  TO  20 

22  WRITE ( 5 , 200) 

200  FORMAT( 2X, 1 9H>0,  GO  ON <0,  END) 

ACCEPT*, Z 
IF ( Z ) 41,41,42 
42  WHITEC5,  250) 

250  F0RMAT(2X,28HNEW  VALUES*  DT,  TMAX,  IPRINT,  F) 

ACCEPT*, DT.TMAX, IPRINT, F 
ISTEP=0 
GO  TO  20 
STOP 
41  END 

FUNCTION  FUNC(TE) 

REAL  NE,MA5E,MASXE, LENGTH, NEP,L 
CC MMON/STOR/XE , PO , C P 1 , CX 1 , NE , NE  P , L , A , R , DT , E , 

+ EP 1 , EX  1 , CEL, MASE,M AS XE,V, LENGTH, PI ,CP1 I ,CPI,EPI 
VD=SQRT( 2*TE/MASXE ) *3E 1 0 
RO6.38E7*(S0RT(2*TE)  )**3 
RP1=RC*EXP(-EP1/TE)*( 1 +EP I /( 2*TE ) )*CPI*PO 
REL= ( CEL*XE*SQRT( 2*TE/MASE ) ) *3E 1 0*( 2*MASE/MASXE ) *TE 
FUNC=( V-R*A*E*NE*VD-L*A*E*VD*(NE-NEP)/OT) /LENGTH 
+ -( 1/VD)*(RPI*EP1+REL) 

RETURN ;END 
FUNCTION  RFN(Y) 

COMMON/STE/TE 

RFN=6 . 3 8E7*( SORT ( 2*TE ) ) **3*E  XP ( -Y/TE )*( l+Y/(2*TE) ) 

RETURN ;END 
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